Market dynamics immediately before and after financial shocks: quantifying the 
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We study the cascading dynamics immediately before and immediately after 219 market shocks. 
We define the time of a market shock Tc to be the time for which the market volatility V{Tc) 
has a peak that exceeds a predetermined threshold. The cascade of high volatility "aftershocks" 
triggered by the "main shock" is quantitatively similar to earthquakes and solar flares, which have 
been described by three empirical laws — the Omori law, the productivity law, and the Bath 
law. We analyze the most traded 531 stocks in U.S. markets during the two-year period 2001- 
2002 at the 1-minute time resolution. We find quantitative relations between the "main shock" 
magnitude M = logj^Q V{Tc) and the parameters quantifying the decay of volatility aftershocks as 
well as the volatility preshocks. We also find that stocks with larger trading activity react more 
strongly and more quickly to market shocks than stocks with smaller trading activity. Our findings 
characterize the typical volatility response conditional on A/, both at the market and the individual 
stock scale. We argue that there is potential utility in these three statistical quantitative relations 
with applications in option pricing and volatility trading. 
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I. INTRODUCTION 

Financial fluctuations have been a topic of study for 
economists P, , mathematicians 0] , and physicists [1- 
[^. Here we study financial fluctuations using concepts 
developed in the field of seismology [13, [HI and analo- 
gies from turbulent dynamics in our description of mar- 
ket main shock magnitudes in order to analyze the dy- 
namic response of markets to financial shocks. We iden- 
tify parallels between energy cascades and information 
cascades, and also between turbulent bursts and the clus- 
tering of volatility [l2|. Our results demonstrate three 
statistical regularities which relate the volatility magni- 
tude M = logiQ ViTc) to the market response before and 
after market shocks. 

Common financial "shocks" are relatively smaller in 
the volatility magnitude, the duration, and the num- 
ber of stocks affected, than the extremely large and in- 
frequent financial crashes. Devastating financial shocks 
such as Black Monday (20 October, 1987) have signif- 
icant aftershocks that can last for several months, and 
this "dynamic relaxation" is similar to the aftershock 
cascade following an earthquake [l^. Here we aim to 
better understand market shocks over a range of M val- 
ues. While the previous studies have focussed on at most 
a few large crashes, we use a large data set of 219 finan- 
cial "main shocks" observed in American markets over 
the 2-year period 2001-2002. We analyze 531 frequently 
traded stocks corresponding to approximately 44,000,000 
volatility records at a 1-minute time resolution. We find 
three quantitative relations which enable answering such 
questions as: 

(i) How does the rate of volatility aftershocks decay 
with time, and how do the decay parameters relate 
to the main shock magnitude Ml 



(ii) How many aftershocks above a given threshold can 
be expected after a main shock of magnitude Ml 

(iii) What is the relation between the value of the main 
shock volatility V{Tc) and the second largest after- 
shock (or preshock)? 

These three questions have been studied for geophysical 
earthquakes, and the corresponding statistical laws are 
referred to respectively as the Omori law, the productiv- 
ity law, and the Bath law. 

The Omori law was first investigated in the context of 
financial crashes by Lillo and Mantcgna [l^ , who found a 
power-law relaxation of fluctuations at a 1-min time res- 
olution for the S&P500 over the 100-day period following 
the Black Monday crash. Power-law relaxation of after- 
shocks is also observed for long periods following several 
other medium-size crashes [l3| , and also for short periods 
up to several days following U.S. Federal Reserve inter- 
est rate change announcements [lH|- One key feature of 
long-range relaxation dynamics is the scale-free decay of 
large fluctuations that is typical of a system with mem- 
ory, and which is complemented by self-similarity in the 
decay substructure |14| . 

We find similar perturbation-response dynamics in the 
intraday volatility (absolute return) time series for many 
single stocks on numerous days, indicating that markets 
respond in a common way to perturbations that range 
in size from everyday market fluctuations to infrequent 
market crashes. Interestingly, the market is very respon- 
sive to Federal Open Market Committee (FOMC) news, 
either in the form of subtle hints from the Fed or ac- 
tual rate changes (expected or unexpected) , because Fed 
Target rates serve as a benchmark and barometer for 
both U.S. and World markets [l^. Methods have been 
developed to use the interplay between the U.S. Trea- 
sury Bill and the Federal Funds effective rate in order 
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to estimate the future movement of the Federal Funds 
target rate [l^. More complex methods to estimate the 
probability of interest-rate change involves anal yzin g the 
price- movement of expiring derivative contracts |17| . The 
connection between macroeconomic factors and financial 
markets is a tribute to the complexity and connectivity 
of economic systems. It is a further indicator that news, 
in addition to complex order-book dynamics, can play a 
significant role in explaining the large rate of occurence 
of large fiuctuations in markets. 

Here we quantify the rate n{\t — Tc\) of volatility shocks 
at time t both before and after a market shock occur- 
ing at time T^. In order to determine Tc, we develop a 
method for selecting a critical time Tc from a set of candi- 
date times {tc} for which the collective market volatility 
of S individual stocks is above a given threshold. For 
19 particular dates corresponding to days with FOMC 
announcements, we compare the values of calculated Tc 
with the reported values of T analyzed in [l^, and we 
find good prediction of T using this method. After this 
calibration, we study the relaxation dynamics of 5 = 531 
stocks, analyzing the Omori law, the productivity law, 
and the Bath law for the dynamics both before (t < Tc) 
and after {t > Tc) the main market shock. 

In Section |ll] we discuss the data, the quantitative 
methods used to calculate n{\t — Tc\), and define col- 
lective market movement. In Section llll Al we quantify 
the threshold for selecting candidate cascades and cali- 
brate using known values of Tc corresponding to FOMC 
meetings. In Section IIII Bl we describe the method for 
choosing Tc from each significant cascade we identify. In 
Section IIVI we discuss the Omori-law parameters a and 
17, the productivity parameter 11, and the Bath law pa- 
rameter B. We note that both 11 and B are indepen- 
dent of the dynamical model, and hence do not depend 
on n{\t — Tel), the functional form of the relaxation dy- 
namics. For each of the statistical laws, we compare the 
results we obtain for the market average with the results 
we obtain for individual stocks. 



II. DATA ANALYZED 

For the two-year period 2001-2002, we analyze Trades 
and Quotes (TAQ) data of more than 500 stocks hsted on 
the NASDAQ and NYSE. In order to analyze the most 
important subset of stocks, we rank each stock by the 
average number of transactions per minute. Wc find 
S' = 531 stocks with an average of more than 3 trans- 
actions per minute, S = 136 stocks with an average of 
more than 10 transactions per minute, and S* = 20 stocks 
with an average of more than 50 transactions per minute. 
Unless otherwise stated, our results correspond to the 
top S = 531 stocks, but all results become more sta- 
tistically significant for smaller subsets of more heavily 
traded (bcUweather) stocks. 

In this paper, we study the volatility Vj(t) of the in- 
traday price time series Pj{t) for stock j. The intraday 



volatility (absolute returns) is expressed as 

v,it)^\Hp,{t)/p,it-5t))\ , (1) 

where here we choose 6t = I minute so that we can an- 
alyze the dynamics immediately before and immediately 
after market shocks. To compare stocks, we scale each 
volatility time series by the standard deviation over the 
entire period analyzed. We then remove the "U" -shaped 
intraday trading pattern (averaged over 531 stocks) from 
each time scries. This establishes a normalized volatility 
in units of standard deviation for all minutes during the 
day and for all stocks analyzed (see Ref. p^). 

We introduce a volatility threshold q which defines a 
binary volatility time series nj(t) for each stock j, which 
we calculate from the normalized volatility time series 
Vj{t) as 

n,it)^(l ' (2) 

We find that a volatility threshold q = 3a is large enough 
to distinguish between significant fiuctuations and nor- 
mal background activity. We also choose this value 
(7 = 3(7 to provide comparison with the analysis per- 
formed in [l^. The rate n{t) measures the fraction of 
the market exceeding q at time t, 

1 ^ 

n{t) = -Y,n,{t) . (3) 

The rate rij (t) quantifying the volatility of a single stock 
j corresponds to the limit 5 — > 1. We define the average 
market volatility V{t) similarly by 

Vit)^^Y.',,{t) . (4) 
* j=i 

A market shock at time Tc may result from exogenous 
(external) news or endogenous herding [isl . [l9| . In many 
cases, the market shocks can be linked to exogenous news 
using archived news feeds that cover and summarize daily 
market events [l^l • In order to analyze market dynamics 
symmetrically around a market shock at time Tc, we an- 
alyze the per unit time rate n(\t — Tc\) around the time 
Tc- It has been empirically observed that the response 
dynamics in financial markets show a power-law decay 
tlSnlSEHll, 

ni\t-Tc\)^a\t~Tc\-'\ (5) 

where is called the Omori power-law exponent, a is the 
cascade amplitude, t < Tc corresponds to before the main 
shock, and t > Tc corresponds to after the main shock. 
For comparison, n{\t — Tc\) is constant for stochastic pro- 
cesses with no memory, corresponding to f2 = 0. Hence, 
for an empirical value w 0, the rate n{\t ~Tc\) is indis- 
tinguishable from an exponential decay for |i — Tcj <t, 
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where t is the characteristic exponential time scale. How- 
ever, for larger values of fl, the exponential and power- 
law response curves are distinguishable, especially if sev- 
eral orders of magnitude in t is available. 

Instead of analyzing n{\t — Tc\) , we perform our quan- 
titative analysis on N{\t — Tc\), the cumulative number 
of events above threshold q at time t minutes, where by 
definition 

N{\t-T,\)^ f n{\t' ~T,\)dt' P{\t~Tc\)^-'' (6) 

for market co-movement and 

NK\t~n)^ f n,{\t' ~T,\)dt' ^ P,{\t-T,\Y'''^ . (7) 

for the activity of stock j. We perform our regression 
analysis on N\\t — Tc\) because it is less noisy and mono- 
tonic as compared to nj{\t — Tc\). 

Hence, for a given day, we calculate the cumulative 
time series N^{t) from nj{t) for each stock j, where t = 
corresponds to the opening bell at 9:30 AM ET. For com- 
parison, we also analyze the average market response 
Nit) of the S stocks under consideration, which com- 
plements the study of individual stocks. 

To demonstrate our approach, in Fig. [T] we plot V(t), 
N{t) and also N^{t) for four single stocks on 01/11/2002, 
a day when there was a large market shock correspond- 
ing to a publicized comment by the Fed chairman Alan 
Greenspan concerning economic recovery which occurred 
at approximately Tc = 255 minutes after the opening 
bell. 

In order to compare the dynamics before and after the 
market shock, we first separate the intraday time series 
N{t) into two time series Nb{t\t < Tc), and Nait\t > Tc). 
Then, to treat the dynamics symmetrically around Tc, we 
define the displaced time t — \t — Tc\ > 1 a.s the temporal 
distance from Tc. As an illustration, we plot in Fig. [2] 
the time scries on 01/11/2002 as a function of r. Wc 
then employ a linear fit to find the r dependence of both 
7Vb(r) = N{Tc) -N{\t- Tc\) and iV,(T) = N{t - Tc) - 
N{Tc) on a log-log scale to estimate the Omori power-law 
exponents fib before the news and fla after the news. By 
analogy, wc define a to be the amplitude a = /3(1 — f2) 
before Tc as af, and after the shock as aa. 

III. METHOD FOR DETERMINING Tc 

A. Calibration using FOMC announcements 

We use n(t) to quantitatively determine times Tc in 
which the market is moving together, possibly in response 
to an external market shock or possibly as a result of en- 
dogenous herding. In Fig. [3] we plot the average daily 
pattern for n{t) and the standard deviation a(t). The val- 
ues of n(t) and (7{t) are not stationary, so we remove the 
daily trend from n{t) by defining the detrended quantity 



n'{t) = ~n{t))/a{t). In order to distinguish signif- 
icant moments of market co-movement from background 
fluctuations, we use a significance threshold which we 
obtain from the distribution of market activity over the 
entire data set analyzed. Hence, we analyze the quantity 
x{t) defined as, 

/ N / N n(t) — n(t)) 

cr(t) 

which is the product of n{t) and n'{t). The value of n{t) 
quantifies the size of the market co-movement, while n'{t) 
quantifies the significance of the market co-movement. 
Because n(t)) is not constant during the day, we con- 
sider the normalized quantity n'{t) in order to remove 
the intraday pattern. Then, to restrict our analysis to 
relatively large market co-movements, we eliminate times 
toward the beginning and end of each day, when average 
market activity is lower (significant morning activity is 
often related to overnight news [l^l). We analyze the 
quantity x{t), which is large only if both n{t) and n'{t) 
are large. Fig. S] demonstrates how the quantity x{t) 
is useful for amplifying market co-movement and pro- 
vides an illustration of a significant shock with substan- 
tial preshock and aftershock dynamics. 

Wc analyze the time series x{t) in order to select the 
set of times {t} of the market shocks that are large in 
the fraction of the market involved (large n{t)) as well 
as significant with respect to the time in which they oc- 
cur (large n'{t)). We determine a significance threshold 
Xc from the probability density function (pdf) of x{t) as 
in Fig. [SJ As a null model, we shuffie the order of each 
intraday time series Vj(t) and obtain a shuffled market 
volatility rate nsh{t) for each day. This preserves the 
empirical pdf of Vj {t) but removes the correlations that 
exist in the temporal structure of Vj{t). We also plot 
nsh{t) « 0.23 in Fig. [3] which corresponds to a residual 
0.23 co-movement due to random fluctuations. We com- 
pare the pdfs for x{t) and Xsh{t) in Fig.[51Jb), and observe 
a significant divergence for x(t) > 1. 

We calibrate our method for determining Tc from can- 
didate cascades by using the known reported values T 
corresponding to Fed announcements. We choose the 
value Xc = 1.0 which reproduces with the best accuracy 
the values of T that we provide for comparison in Ta- 
ble H] The value of Xc = 1.0 results in 5,804 minutes 
out of 190, 000 minutes analyzed for which x{t) > Xc, or 
roughly 3% of the 2-year period with significant market 
CO- movement. 



B. Determining Tc from candidate cascades 

In a typical trading day there are many large fluctua- 
tions, for both individual stocks and indices such as the 
S&P 500 and DOW. Analysis of the return intervals and 
the cross-correlations in financial time series shows that 
significant statistical regularities exist [25-31]. This fact 
is evident in the robust probability density function of 
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TABLE I: Comparison of announcement times T (as re- 
ported in New York Times) with the market clustering times 
Tc calculated using a threshold Xc ~ 1.0, cascade window 
Al = 60 min., and S = 136 stocks. The value of x{Tc) cor- 
responds to the largest value out of all the candidate {x} in 
the most significant cascade of the particular day. Dates of 
19 FOMC meetings in the 2-year period between Jan. 2001 
- Dec. 2002, where the Federal Funds Target rate (Rnew) 
was implemented by the rate change (AR) at (T) minutes 
after the opening beU at 9:30 AM ET [32]. The absolute rel- 
ative change |-^| = \AR{t)/R{t- 1)| has typically filled 
the range between 0.0 and 0.25. Note: Date** refers to un- 
scheduled meetings, in which the announcement time did not 
correspond to 2:15 PM ET (T = 285 minutes). 



^OMC Date 


Rnew (%) 


AR 


AR 


T 


Tc 


01/03/01** 


6 


-0.5 


-0.077 


210 


227 


01/31/01 


5.5 


-0.5 


-0.083 


285 


290 


03/20/01 


5 


-0.5 


-0.091 


285 


286 


04/18/01** 


4.5 


-0.5 


-0.100 


90 


88 


05/15/01 


4 


-0.5 


-0.111 


285 


287 


06/27/01 


3.75 


-0.25 


-0.063 


285 


285 


08/21/01 


3.5 


-0.25 


-0.067 


285 


286 


09/17/01** 


3 


-0.5 


-0.143 





16 


10/02/01 


2.5 


-0.5 


-0.167 


285 


288 


11/06/01 


2 


-0.5 


-0.200 


285 


292 


12/11/01 


1.75 


-0.25 


-0.125 


285 


287 


01/30/02 


1.75 





0.00 


285 


289 


03/19/02 


1.75 





0.00 


285 


293 


05/07/02 


1.75 





0.00 


285 


287 


06/26/02 


1.75 





0.00 


285 


286 


08/13/02 


1.75 





0.00 


285 


291 


09/24/02 


1.75 





0.00 


285 


291 


11/06/02 


1.25 


-0.5 


-0.286 


285 


286 


12/10/02 


1.25 





0.00 


285 


295 



volatility which has a stable power-law tail for a wide 
range of time scales ranging from 1-minute to several days 
d, 0, 113, nil . We select market cascades that are above 
a "spurious fluctuations" threshold, which we define by 
randomizing the order Vi{t). We use the corresponding 
shuffled values Ush (t) as a proxy for background noise. 

We find on average approximately 12 minutes per day 
above the threshold Xc = 1.0. So here we develop a 
method for selecting the most likely time Tc from all can- 
didate times with x{t) > x^. For a given day, we collect 
all values of x(fy > Xc into a subset {x'{t)} of size z. 
From this subset, we divide the z values into k cascades 
{x'{t)}i, which we define as localized groups using the 
criterion that a cascade ends when the time between the 
last x' in cascade i is separated from the first x' in cas- 
cade i + 1 by a time window greater than Al minutes. We 
next assign to each cascade group {x'{t)}i a weight equal 



to the sum of the x'{t) values belonging to the given cas- 
cade group, and select the cascade group with the largest 
weight as the most significant cascade. Within the most 
significant cascade group, we choose the time correspond- 
ing to the maximum value of x' (t) as the time Tc of the 
main shock. We calibrate this method using the reported 
times for the 19 FOMC interest rate meeting announce- 
ments, and find that the values Al = 60 and Xc = 1.0 
best reproduce the known set {T}, which we provide for 
comparison in Table HI 

Using the parameter values Xc ~ 1.0 and Al = 60, 
we find 373 days with market shocks, out of 495 days 
studied. If the values of x'{t) were distributed uniformly 
across all days, then the probability of finding 122 days 
without one x'{t) is vanishingly small, which confirms 
that the x'{t) group together forming cascades. We re- 
move all days where Tc is within At = 90 minutes of 
opening {t = 0) or closing {t = 390), and all Tc that 
occur on half-days (days before or after the 4th of July, 
Thanksgiving or Christmas), resulting in the data set 
{Tc} constituting 219 individual days. We also analyzed 
the subset of 156 market shocks within At = 120 minutes 
of opening or closing and find analogous results as those 
reported here for At = 90. 

Furthermore, in order to test the dependence of the 
data set {Tc^^} found for the time resolution St = I 
minute used in this paper, we also compare the values 
of{ri^^}and {Ti'°'} found using a volatility series with 
dt = 5 min. and 6t = 10 min. resolution, respectively 
(sec Eq. ([T|)). For each of the 219 days with a Tc value we 

calculate the absolute difference in the time value Tc^*^ 
using two values of St. We use similar values of Xc for each 
time resolution so that the number of days with market 
shocks for each resolution are approximately equal. The 
difference in ri'^*'' depends on the resolution St and the 
locality STc associated with each market shock. The av- 
erage of the absolute differences for three values of St are 
_ j.a)| ^ g nrinutcs and |Ti^"^ -Ti^'| = 15 minutes. 
We estimate the standard error for a particular time reso- 
lution STc^^^ « 2St, which implies that STc^^ « 2 min. for 
the 1-min. time resolution. Hence, the use of Tc ± STc 
does not significantly change the results of this paper. 
In the next Section, we analyze the empirical laws that 
quantify the response dynamics both before and after 
significant market shocks. 



IV. RESULTS 

The analysis performed in this paper is largely inspired 
by the analogies between financial market crashes and 
earthquakes. Here we identify 219 cascades that meet 
our volatility significance criterion and that are within 
At = 90 minutes of the open or closing of the trading 
day. For this set of shocks, we analyze the dynamics 
over the 90 minute period immediately before Tc and 
the 90 minute period immediately after Tc using the 
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framework developed in earthquake research [39|-|45| . Al- 
though we present results for only At = 90, we also an- 
alyzed the subset of 156 shocks within At — 120 min- 
utes of the beginning and ending of the trading day, and 
find analogous results. We restrict our analysis to the 
local time period 2At within the trading day so that 
our results arc minimally affected by overnight effects 
and overlapping shocks (since the frequency of shocks 
larger than our threshold Xc used here is approximately 
one per trading day). Hence, our analysis of the size- 
dependence of market relaxation dynamics, where we re- 
late the cascade parameters to the market shock magni- 
tude M = \ogiQ V{Tc), pertains to the intraday horizon 
of market shocks which occur relatively frequently. 

A recent study finds significant evidence of Omori 
power-law relaxation both before and after common 
FOMC interest rate announcements [l5|, and uses the 
relationship between the overnight Effective rate and the 
U.S. 6-month Treasury Bill to estimate the magnitude 
of the financial news shock. The dynamics before the 
Fed announcements, which are regularly scheduled and 
pre- announced, are consistent with market's anticipated 
surprise in the Fed news, while the dynamics after the 
announcements are related to the market's perceived sur- 
prise in the Fed news. The Federal Reserve uses the "an- 
nouncement effect" [i^ to control the federal funds mar- 
ket. Despite the calculated monetary policy, the markets 
react quite strongly to the news, with large Omori-law 
relaxation cascades that typically correspond to the rel- 
atively large M found in this paper. 

Closely related to the Omori relaxation of aftershocks 
is the productivity law, which establishes a power- 
law relationship between the number of aftershocks or 
preshocks that follow or precede a main shock of mag- 
nitude M. This is analogous to earthquake analysis 
where the magnitude is defined as M « (2/3)logiQi? 
j43| . where E is the energy associated with the stress 
released by the main shock. We justify our analogy be- 
tween market volatility V and earthquake energy E by 
comparing the cumulative distribution 

P{V > s) - s-"^ (9) 

of volatility in financial markets with the cumulative dis- 
tribution 

P{E > s) - s-''^ (10) 

of energy E in seismic earthquakes. Both cumulative 
distributions are asymptotically power laws, with ny w 3 
milil and the Gutenberg-Richter law 7?^ « 2/3 [ll. 

For the set of 219 market shocks we analyze, we find 
a wide range of V{Tc), and hence a wide range of cas- 
cade dynamics. We analyze the dynamics only within the 
trading day so that we can be confident that the dynam- 
ics after and before Tc are related to the market shock 
V{Tc). Also, for relatively small V{Tc), it may be diffi- 
cult to distinguish cascade preshocks (aftershocks) from 
opening and closing effects, and volatility resulting from 



overnight news. Furthermore, we only analyze the first 
At = 90 minutes of each Na{T) and iVb(r) time series, 
so that the comparison of productivity Pa_f,(At) is not 
affected by time series of variable length. 

In Fig. [S] we plot the pdf of Omori parameter values 
^a,b and aa,b obtained from the power-law fits of A^6(t) 
and Nair). Figs. and|6l[b) show the distribution of 
parameter values calculated for the average market re- 
sponses Nh{T) and Nair) corresponding to Eq. while 
and ^d) show the distribution of parameter values 
calculated from the individual stock responses Nl (r) and 
N-l^ (r) . The pdfs for individual stock values of V, and a 
have a larger dispersion, as the response to each market 
shock is not uniform across all stocks. For the average 
market response N{t) in Figs. [6ja) and[6jb) the pdfs of 
il and a are shifted to larger values for t > Tc as com- 
pared to t < Tc- This is indicative of the stress that can 
build prior to anticipated announcements (46j and the 
surprise that is inherent in the news. Larger n values 
correspond to faster relaxation times, while larger a val- 
ues correspond to higher activity. We also observe £7 < 0, 
which corresponds to particular time series in which the 
pre-shocks or after-shocks farther away from the main 
shock (for large r) are dominant over the volatility cas- 
cade around Tc- The values of the Omori parameters we 
find on averageing over all market shocks arc given in the 
figure caption. 

Although there is a wide distribution of Omori param- 
eter values when considering all 219 market shocks, there 
is a strong correlation between the individual stock dy- 
namics for a given market shock. In Fig. [7] we relate 
the values of a and calculated for the average mar- 
ket response to the average and standard deviation of 
a and 51 calculated for individual stocks for a given Tc- 
The strong correlation between these quantities over 219 
different dates indicates that the dispersion in the val- 
ues of a and for individual stocks, as demonstrated in 
Figs, inta) andin^c), result from the broad range of V{Tc) 
values, and further, that the dispersion does not result 
merely from the range of stocks analyzed. 

In Fig. [5] we plot the relation between the magnitude 
M of each main shock and the resulting Omori exponents 
51a. f, calculated from both market Na,b{T) and individ- 
ual stock jj(t) response curves. Figs. [Slja) and [HJc) 
show a positive relation between M and the decay expo- 
nent Qa, which indicates that the market responds faster 
to large shocks on the intraday time scale. Figs, ^h) 
and [HKd) show a significant dispersion across all stocks 
for a given date. Interestingly, we find a crossover at 
A/x ~ 0.5 above which Qa.t increases sharply to positive 
values. The values 51 w for M < Mx correspond to 
a dynamical cascade ^(t) that is indistinguishable from 
an exponential decay. Typically, small values of 51 corre- 
spond to stocks with relatively low trading activity which 
are less sensitive to market shocks. For individual stocks, 
we define M to be the logarithm of the largest volatil- 
ity within 3 minutes of the main shock Tc measured for 
the average market response N{t). This accounts for the 
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possibility of a stock-specific anticipation or delay time in 
the volatility as a result of the mainshock V{Tc). There 
is also the possibility that a spurious value of 51 w 
can arise from a stock which has a constant (potentially 
large) level of activity throughout the entire time period 
analyzed. 

In Fig. [9] we plot the relation between the magnitude 
M and the Omori-law amplitude aa.t for both market 
Na,b(T) and individual jj(r) response curves. The re- 
lation between a and M is stronger than the relation 
between f2 and M, indicating that the Omori-law ampli- 
tude has a higher information content than the Omori- 
law exponent. The strong relation for the average market 
response suggests that it is possible to identify precur- 
sors of market shocks with statistical certainty. How- 
ever, since often T^. corresponds to anticipated market 
news, the significant activity prior to the main shock is 
a natural biproduct of trader anticipation. Interestingly, 
we also observe a critical threshold for Mx « 0.5, above 
which the average response amplitude aa.b increases sud- 
denly, analogous to a first order transition. 

The behavior of the market cascades above and below 
A/x are significantly different. Below Mx, it is typical 
for fla.b to be negative and aa,b ~ 0, indicating that 
the triggering shock V{Tc) is relatively insignificant, with 
relatively few preshocks and aftershocks. In this scenario, 
it is possible for a few clusters of relatively large volatility 
towards the end of the time series ^^(t) to dominate 
in the calculation of the best-fit parameters, producing 
negative values for fla.b- There are also cases for both 
j(t) and Na^ir) for which ila.b ~ 0, corresponding to 
a constant rate of preshocks and aftershocks in the time 
period analyzed. 

However, above Mx, the cascade around Tc is more sig- 
nificant, but with some anomalous properties. Namely, in 
the case of the the market response Na,biT), there is an in- 
creasing relation between AI and Cla , indicating that the 
market responds more quickly to larger market shocks. 
This observation is consistent with the "semi-strong" 
efficient-market hypothesis, which asserts that markets 
incorporate public news approximately instantly. This 
observation is similar to geophysical earthquakes, where 
it is observed that the value of Via increases with M for 
a given geographical location [13]. Also, for geophysical 
earthquakes, it is also found that the Omori amplitude 
increases exponentially with M [42| . In the case of indi- 
vidual company response NI^^{t), the values of the av- 
erage ^la^b ~ 0.1 saturate above Mx, whereas the values 
of aa.b increase with M. Thus, the stochastic dynamics 
display a non-linear relationship with A/, consistent with 
a non-linear shot noise process [48| . and are a potential 
avenue for future theoretical research. 

In Fig. [TU] we plot the relation between V{Tc) and 
the productivity Pa{At) (or Pb{At)), defined as the cu- 
mulative number of aftershocks (or preshocks) greater 
than the threshold q = 3 within At = 90 minutes of 
Tc- Motivated by the power-law relationship observed 
for earthquakes we fit the relations Pa{At) ^ M^'^ and 



Pb{At) ^ M^'', and find statistically significant values 
for the market response H;, = 0.38 ± 0.07 and Ha = 
0.48±0.04, and for individual stocks H,, = 0.23±0.01 and 
Ha = 0.25±0.01. For earthquakes, [i^ reports a range 
of Ha ?» 0.7 — 0.9 values that are larger than observed 
here for financial markets, meaning that the productiv- 
ity of physical earthquakes increases "faster" with main 
shock magnitude than does the productivity of market 
shocks. Since for earthquakes Ha < &, where & « 1 is the 
scaling exponent in the Gutenberg-Richter law, this in- 
equality establishes the relative importance of small fluc- 
tuations as compared to large fluctuations [l^] . In other 
words, this inequality indicates that small earthquakes 
play a larger role than large earthquakes in producing 
the observed number of large earthquake shocks. Using 
an analogous argument for market volatility, since the 
cumulative distribution exponeiit riy « 3 is found to be 
robust across many markets [stI [3a | , then the total num- 
ber NTot{V) of aftershocks triggered by all main shock 
of size V scales as, 

NTotiV) = P{V)Pa{At) ^ lo(n»-w)iogv ^ (^^^ 

is a decreasing function of V. Hence, we also flnd that 
aftershock cascade triggering is controlled by the con- 
tributions of the more numerous small shocks. Thus, 
the medium-sized market shocks (analyzed here) play a 
larger role than the large market shocks in producing the 
observed heavy-tailed distribution of market shocks. We 
further note that the productivity is a combination of the 
relationships of both a and with V{Tc), which can be 
written as 

Pa{At) = Na{At) ^ {At)^-''- aa/{l - ila) ~ V{T,f^ , 

(12) 

with equivalent relation before the shock for Pf,(At). 

In Fig. [11] we plot the values il, a, and P{At), both 
before and after the main shock at time Tc- Surprisingly, 
while there is little statistical relation between fif, and 
51q, there is a strong relation between ah and aa as well 
as between Pb{At) and Pa{At), for both At 90 and 
At = 120 minutes. This result could be of interest for 
volatility traders and options traders who would like to 
anticipate the market dynamics after an announcement, 
given the dynamics before the announcement. 

In Fig. [12] we relate the size of the largest shock Vi = 
V{Tc) to the sizes of the second largest shock V2, both 
before and after Tc- The Bath law parameter B quantifles 
the relation between Vi and V2 as 

Ml - Ah = log Vi - log V2=B . (13) 

This functional form implies the relation 

V2IV1 = Cb (14) 

and hence B = — log^o Cb- Fig. W2\c) is a scatter plot of 
Vi and V2,a which shows a linear relation corresponding 
to Ba = -^logio(0.90) = 0.046. Surprisingly, Fig. [HJa) 
also shows a strong relation between Vi and V2,6 with 
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Bb = — log]^g(0.81) = 0.092. Comparing the values of 
and Sq, the difference between the Vi and V2 is smaller 
after Tc than before Tc- Interestingly, both Bi, and Ba 
are significantly smaller than the value Be ~ 1.2 ob- 
served for earthquake aftershocks [4l|, meaning that the 
largest preshock and aftershock are of comparable mag- 
nitude to the main shock. This significant difference be- 
tween earthquakes and market shocks is largely due to 
the relative probabilities of observing first and second- 
largest events xi and X2- The conditional probability 
P(xi\x2) = P(xi > X2) is given by the corresponding 
cumulative distribution function. Hence, using Eq. Q 
and Eq. ^TU\\ . the ratio of the conditional probabilities 
for El and Vi is 



PiVl\V2. 

P{Ei\E2 



PjVi > V2) 
P{Ei > E2) 



E-^l^ 



(15) 



which roughly explains the 10^ factor difference Be ~ 
WBv 

We also compare the volatilities V\ and V2 for individ- 
ual stocks in Fig. [T^b) before and in Fig. [T^d) after 
Tc. Wc compute the average value (V2) for linear bins, 
and find V\ > {V2) for Vi > 20, both before and after Tc. 
Also, Fig. [12] shows that (1^2,0) > {V2,b) for most values 
of V{Tc). Hence, the reaction to surpise causes larger 
volatility fluctuations than the anticipation of surprise. 

We further ask the question, how do the response pa- 
rameters analyzed here depend on the variations between 
individual stock trading patterns? To answer this ques- 
tion, we quantify the trading capacity of each stock by 
(w), the average number of transactions per minute, with 
3 < (w) < 163 for the S = 531 stocks analyzed. Wc hy- 
pothesize that (w) is closely related to firm size and mar- 
ket impact. FigfT3la) shows that (a), (w) and {P{At)) 
increase with (lj) after Tc, indicating that stocks with a 
large trading base respond to market shocks with large 
volatility {v{Tc)) (shown in Fig lTST b)). but also relax 
more quickly, corresponding to larger D, values. However, 
we find no statistical relation between (w) and {v2,a)- In- 
terestingly, Fig. [TSf b) shows that this positive relation 
also applies to the dynamic response parameters before 



V. DISCUSSION 

Cascading avalanche dynamics are a common phenom- 
ena in complex systems rangiiig in scale from solar flares 
[39I lioj and earthquakes fiol [Til lisl . to microscopic 
vortices in turbulent fluids |49j . Similar bursting phe- 
nomena is also observed in human organs, such as the 
heart [10, [5l|, lun gs Isp. 531. and brain (s^^SGj . and also 
for common social [1^57-53] and economic systems [12I - 
[Tsl . [2l| - |23l . [60I HH. Neural avalanches in the brain are 
frequent even in the resting state, and are a signature of 
healthy brain functioning within the neural network. In 
fact, the ability to process and disseminate information 



is largely attributed to the network structure of neuronal 
correlations which, if inhibited by disease, lead to altered 
disfunctional states such as in the case of schizophrenia. 
Extending by analogy, the frequency of cascades in fi- 
nancial markets could also be viewed as a "healthy" op- 
timal state for processing information and eliminating 
arbitrage among the many the degrees of freedom. Re- 
cent work [62j on the switching dynamics around highs 
and lows in finanancial time series shows further evidence 
of Omori power-law scaling before and after microtrend 
extrema, in analogy to the market shocks at Tc developed 
here. Interestingly, this work on switching dynamics finds 
cascading trends on time scales ranging from seconds to 
hundreds of days. 

Financial markets are subject to constant information 
fiow, resulting in a lar ge r ate of significant events, such 
as Fed announcements [isl |46| . quarterly earnings, splits 
and dividends announcements, mergers and acquisitions, 
institutional reports. This information can arrive as "ex- 
pected" or come as a "surprise" . Interestingly, there are 
precursors extending more than a day in advance of ex- 
pected announcements such as earnings announcements 
[6(j |. Economists have long been interested in the inter- 
play between informed and uninformed traders, and the 
dissemination of information across a market consisting 
of rational agents. Early work focusses on the relation- 
ship between trading volume and price change, and the 
relationship between these quantities and the qualitative 
notions of surprise, importance, and precision of the in- 
formation [6l| . 

Using methods from statistical physics [4,5,34] and 
geophysics, we analyze the absolute returns of price be- 
cause of the long-memory property, and the universality 
of fluctuations in this quantity across diverse markets 
[35I . [36l [ssj . Ref.[6l| postulates that price changes re- 
flect the average change in market expectations, whereas 
trading volume reflects idiosyncratic reactions across all 
traders. Recent work further quantifies trading volume 
fiuctuations and finds that they are similar to price 
fiuctuations, and furthermore, finds significant cross- 
correlation between volume change and price change j63| . 
Omori relaxation dynamics are also shown for trading 
volume in [isj . Here we also observe significant volume 
cascading as evident in Fig. 2] The analysis of volume 
and transaction dynamics is an avenue of future research, 
and could highlight the relationship between volume and 
price fiuctuations by studying their correlation around 
market shocks. 

To summarize, we analyzed the cascade dynamics of 
price volatility, which has potential applications in op- 
tions pricing and the pricing of other derivatives. The 
Black-Scholes equation in its simple form assumes that 
the fluctuations in the price are constant during the du- 
ration of the option [6^. However, more sophisticated 
methods 65[ incorporate time dependent price volatility, 
and are more realistic descriptions of the non-stationarity 
of financial time series. The results in this paper are 
of potential interest for traders modeling derivatives on 
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short time scales around expected market shocks, e.g 
earnings reports. The statistical regularity of both mar- 
ket and individual stock behavior before and after a mar- 
ket shock of magnitude M = logj^Q V{Tc) provides in- 
formation that could be used in hedging, since we ob- 
serve a crossover in the cascade dynamics for M « 0.5. 
Knowledge of the Omori response dynamics provides a 
time window over which aftershocks can be expected. 
Similarly, the productivity law provides a more quan- 
titative value for the number of aftershocks to expect. 
Finally, the Bath law provides conditional expectation 
of the largest aftershock and even the largest preshock, 
given the size of the main shock. Of particular impor- 
tance, from the inequality of the productivity law scaling 



exponents and the pdf scaling exponent for price volatil- 
ity, we find that the role of small fluctuations is larger 
than the role of extremely large fluctuations in account- 
ing for the prevalence of aftershocks. 
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FIG. 1: (color online) Typical volatility curves on 01/11/2002 with market shock at Tc = 256 minutes, (a) The cumulative 
volatility N-' (t) for the stock of several large companies have varying behavior before Tc, but each stock shown begins to cascade 
soon after Tc- The market average N{t) over all S — 531 stocks analyzed demonstrates a distinct change in curvature a.t t — Tc- 
(b) The average market volatility V{t) demonstrates a sharp peak at Tc, and also two precursor events at t « 190 and « 230 
minutes. 
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FIG. 2: (color online) (a) An illustration of Ni,{t) and Nair) for the same set of curves plotted in Fig. [T] The displaced time 
T — \t ~ Tc\ is defined symmetrically around Tc = 256 minutes on 01/11/2002. (b) logj^Q A''6(r) and log-^Q Nair) are linear 
with logj^o T over two orders of magnitude on a logarithmic scale. The Omori parameters in Eq. ([5]) calculated from N{t) are 
fli = 0.09 ± 0.01, at = 0.21 ± 0.01 and Qa = 0.32 ± 0.01, Oa = 0.81 ± 0.01 . 
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FIG. 3: (color online) The fraction n{t) of the market above the volatility threshold q is non-stationary through the trading 
day. We plot in (a) the average daily trading pattern n{t) for S = 136 stocks and in (b) the corresponding standard deviation, 
to demonstrate the trends we remove in the normalized quantity n'{t). In practice, we use the smoothed average of these 
curves in order to diminish statistical fluctuations on the minute-to-minute scale. For comparison, we compute nsh(t) ~ 0.23 
and ash ~ 0.09 for shuffled Vi{t). The values of n{t) provide an estimate for the background market co-movement that can be 
attributed to random fluctuations. 
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FIG. 4: (color online) Example of market co-movement n{t) in both price volatility and total volume, and the qualitative 
relationship between the quantities n{t), n'{t) = {n{t) — n{t))/a{t), and x(t) = n{t)n'{t). The market shock on 01/11/2002 
occurred at Tc = 256 in response to a public comment by the Fed chairman A. Greenspan concerning economic recovery [20l |. 
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FIG. 5: (color online) Using the volatility threshold g = 3 and S — 136 stocks, we determine the market comovement threshold 
Xc from the pdf of x{t) = n{t)n'{t). (a) The pdf for the 190,000 minutes analyzed of the volatility rate n{t) corresponding 
to the fraction of the market with volatility Vi{t) > q. (b) The pdf of x{t), where in the quantity x(t) we have removed the 
average daily trend of n(t), so that x(t) is relatively large when market comovement is large and significant. For comparison, 
we also plot the pdf of Xghit) computed from randomly shuffled volatility time series Vi{t). We find a divergence between the 
pdf of x{t) and of Xsh(t) for x > 1.0, which we define as the comovement threshold a::c = 1 in our analysis. 
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FIG. 6: (color online) (a,b) Comparison of the probability density functionss P{a) and P(S1) of Omori parameters a and Vl 
computed from the average market response Na,b{T). (c,d) The analogous pdf plots computed from individual stock response 
bi''')- The average and standard deviation of each data set are (a) fla ~ 0.09 ± 0.07 , fib = 0.06 ± 0.07 (b) aa ~ 0.35 ± 0.11 
, a'b = 0.28 ± 0.09 (c) ila = 0.08 ± 0.20 , fib = 0.03 ± 0.22 and (d) = 0.53 ± 0.25 , at = 0.46 ± 0.24. Values of both Qa and 
aa are consistently larger than Qb and at, indicating that the response time after Tc is shorter than the activation time leading 
into Tc. However the response cascade after Tc has, generally, a larger amplitude. 
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FIG. 7: (color online) In order to account for the dispersion in the pdfs plotted in Figs. (Sfc) and|6jd) for individual stocks, we 
compare the average values aa,b and ila,b computed from all A'';^ ^(r) with the aa,b and Q,a,b computed from the corresponding 
average market response Na,b{T) for each of the 219 Tc. The visually apparent correlation indicates that the parameters 
quantifying Na,b{T) are a good representation of the average ^(t). The correlation coefficient r for each linear regression is 
provided in each panel. 
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FIG. 8; (color online) The relation between the magnitude M = logjQ V{Tc) and the Omori exponents i}a,b- In panels (a) 
and (c) we compare values calculated from the average market response Na,b{T), and in panels (b) and (d) we compare values 
calculated from individual stock response A'^;^ (,(''")■ (a-) Weak relation before Tc, where we validate the linear regression model 
at p = 0.001 significance level, but with correlation coefficient r — 0.22. The dispersion may result from the variability 
in anticipation preceding the market shock at Tc. (c) The relation between ila and M is stronger after Tc than before Tc, 
with linear regression significance p ~ 0, correlation r = 0.40, and regression slope m = 0.19 ± 0.03. The increasing trend 
demonstrates that a faster response, quantified by larger fia, follows a larger M. Data points in panels (a) and (c) denoted by 
the symbol x correspond to values of Sl^.b calculated for randomly selected Tc on those 118 days analyzed without a single value 
of x{t) > Xc- In panels (b) and (d) there is much dispersion in the SI values of individual stocks for given V{Tc). However, the 
average trends demonstrate a significant crossover at Mx « 0.5 from ^la,b < to Qa.b > 0. The case of f2 < can occur when 
there is more volatility clustering for large r than for small r, whereas the case of f2 > occurs for large volatility cascading 
around r > 0. This crossover could result from the difference between anticipated and surprise shocks at Tc- 
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FIG. 9: (color online) The relation between the magnitude M = logj^Q V{Tc) and the Omori amplitudes Qa,(,. In panels (a) 
and (c) we compare the values calculated from the average market response N{t) and in panels (b) and (d) we compare values 
calculated from individual stock response (,(''")• (a) The increasing relation between at and M is statistically stronger 
than the relation between Jli, and M in Fig. EJa), with significance p « 0, correlation coefficient r — 0.52 and regression 
slope m — 0.35 ± 0.04. (c) The relation between Oa and M is strong, with significance p ~ 0, r — 0.84, and regression slope 
m = 0.68 ±0.03. Data points in panels (a) and (c) denoted by the symbol x correspond to values of aa,b calculated for randomly 
selected Tc on those 118 days analyzed without a single value of x{t) > Xc- The result that a increases with increasing V{Tc) 
holds even for random times. In panels (b) and (d) there is much dispersion in the a values of individual stocks for given 
V{Tc). However the average trends demonstrate a significant crossover at Mx « 0.5 from aa.b ~ 0.2 for M < 0.5 to Oa.b > 0.2 
for M > 0.5. This crossover occurs at a similar location as the crossover observed in Figs. [Sjb) and (d) for Q.b,a- The average 
amplitude value a increases sharply for M > Afx, consistent with first order phase transition behavior. 
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FIG. 10: (color online) The increasing relation between the productivity Pa.b(At) of each market shock and the size of the 
main shock A'l = logj^g ^(Tc) with At = 90 min. As is found in earthquakes, we find a power-law relationship between M 
and V{Tc) described by a productivity exponent 11;, before and exponent Ha after the market shock. Data points in panels 
(a) and (c) denoted by the symbol x correspond to values of Pa^biAt) calculated for randomly selected Tc on those 118 days 
analyzed without a single value of x{t) > Xc- The result that P{At) increases with increasing V{Tc) holds even for random 
times. For the average market response Nb,a{At), we find (a) lit = 0.38 ± 0.07 and (c) Ha = 0.48 ± 0.04. For the productivity 
of individual stocks corresponding to ^{At) we find (b) 11^ — 0.23 ± 0.01 and (d) Ha = 0.25 ± 0.01. For comparison, the 
power-law exponent value pertaining to earthquake aftershocks is Ua ~ 0.7 — 0.9 [43 |. 
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FIG. 11: (color online) A comparison of Omori parameters before and after Tc for N{t) and varying At indicate that Q(, and 
Pi,{At) axe better conditional estimators for the dynamics after Tc. (a) Weak relationship between f2(, and Qa for At — 90 
and 120. (b) Strong relationship between qj, and Oa for At — 90 and 120, with both linear regressions passing the ANOVA 
F-test at the p < 0.001 confidence level, (c) Strong relationship between P(,(At) and Pa{At) for At = 90 and 120 min. at the 
p < 0.001 confidence level. 
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FIG. 12: (color online) The increasing relation between the size of the main shock V(Tc) and the size of the second largest 
aftershock (or preshock) V2(Af) within At minutes of Tc demonstrates that the volatility of the largest aftershock (or preshock) 
increases with mainshock volatility. As with the Bath law for earthquakes, we observe a proportional relation V2,a(Af) = 
CbV{Tc) which corresponds to a Bath parameter B = — logjgCB. For the average market response Nb,a{A.t) we calculate 
Cb for (a) the dynamics before, Cb ~ 0.81 with correlation coefficient r = 0.70 and ~ 212, and for (c) the dynamics after 
Cb ~ 0.9 with r = 0.87 and ~ 109. For the Bath law corresponding to individual stocks we find that a linear function best 
fits the relation between V{Tc) and the average value V2{At) calculated for equal-sized bins as indicated by circles with one 
standard deviation error bars. We calculate the regression slope for the Bath law (b) before is m = 0.65 ± 0.02 and (d) after is 
m = 0.40 ±0.01 




FIG. 13: (color online) Relations between individual stock trading activity and dynamic response parameters (a-d) after Tc and 
(e-h) before Tc, averaged over all the days with a market shock. We measure the trading activity (uj) for each stock, defined 
as the average number of transactions per minute over the 2-year period 2001-2002. We find that stocks with large trading 
activity react both more strongly (larger a and larger P{At)) and quickly (larger SI) to market shocks. However, panel (d) 
shows that there is little relation between (a;) and the average size of the largest aftershock («2). 



